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Introduction 


The  potential  flow  around  asymmetrical  bodies 

in  the  presence  of  a  free  surface  is  of  interest  in 

many  problems  of  ship  hydrodynamics.  Asymmetry  in 

ships  can  either  be  permanent  as  in  catamarans,  or 

temporary  due  to  rotations  of  symmetric  hulls  in 

yaw  or  roll.  The  purpose  of  this  work  is  to  investigate 

theoretically  the  forces  and  moments  acting  on  such 

bodies  for  the  case  of  a  steady  translational  motion 

in  an  ideal  fluid.  The  solutions  are  limited  to  thin 

bodies  with  a  small  degree  of  asymmetry.  These  restrictions 

are  employed  so  as  to  yield  a  well  defined  boundary-value 

problem  having  a  solution  which  lies  within  the  scope 

of  the  potential  theory. 

The  treatment  of  hydrodynamic  problem  was  based 

mainly  on  information  deduced  from  the  theory  of  low- 

aspect  ratio  wings.  Davidson  and  Schiff  (1946)  pointed 

out  that  larger  changes  in  the  wave-making  pattern  are 

observed  when  the  Froude  number  is  larger  than  C.19. 

Hu  (1961)  solved  this  problem  for  a  yawed  ship  by  the  use 

of  an  iterative  perturbation  method  based  on  the  asymptotic 

expansion  of  the  Kernel  function  of  the  integral  equation 

for  small  Froude  numbers.  He  found  that  the  magnitudes 

of  the  forces  and  moments  acting  on  the  ship  increase 

rapidly  as  the  Froude  number  increases  up  to  F  *  0.35  and 

n 

then  remain  mainly  constant. 


The  solution  of  the  boundary-value  problem  leads 
to  a  Fredholm  Integral  equation  of  the  first  kind  for 
the  doublet  moment  distribution  similar  to  the  usual 
representation  of  a  lifting  surface.  Due  to  the 
complicated  form  of  the  Xernel  function  approximate 
numerical  methods  must  be  applied  to  solve  this  integral 
equation. 

In  this  work  a  numerical  solution  is  used  based 
on  the  finite-element  method,  where  the  region  of 
integration  is  divided  into  small  rectangular  elements 
and  the  unknown  doublet  moments  are  defined  by  an 
approximating  function  within  each  element.  In  doing 
that  and  performing  the  integrations  within  each  element 
analytically,  the  integral  equation  reduces  to  a  set 
of  linear  algebraic  equations  that  can  be  solved  for  the 
unknown  doublet  strength.  We  carry  through  this 
procedure  for  two  cases,  one  when  l/'J?" 20  and  one  when 
L/T-7.  The  results  are  summarized  in  the  concluding 
section  of  this  work. 


I .  Mathematical  Formulations 


It  will  be  convenient  in  formulating  this  problem 
to  introduce  three  right-handed  coordinate  systems. 

One  is  fixed  in  space,  Oxyz,  with  Oy  directed  oppo¬ 
sitely  to  the  force  of  gravity,  Ox  coincides  with  the 
direction  of  motion,  and  Oxz  lies  in  the  plane  of  the 
undisturbed  free  surface.  Of  the  other  two  coordinate 

t  t  t  t 

systems,  one  is  fixed  in  the  body,  Oxyz,  and  one  is 

moving  with  the  body,  Oxyz,  but  is  taken  in  such  a  way 

that  the  (x,z) -plane  coincides  with  the  mean  water 

surface,  and  Ox  makes  an  angle  6  with  Ox.  Further, 

1 1 1 1 

Oxyz  and  Oxyz  coincide  when  the  body  is  at  rest  (see 
Figure  1) . 


x ' 


Figure  1 

Plan  View  of  an  Asymmetric  and  Yawed  Body 
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The  hull  form  will  be  described  by  the  surface 
or  more  explicitly  by 

a'  »  /  (x',y'i  for  the  starboard  surface 

I  0 

and  a'  -  fp(x',y')  for  the  port  surface. 

Define  f  (x'ty')  and  ft(x',y')  by 
-  1/2  U  (x\y')  +  f  (*',V ') ) 
ft(x',y')  -  1/2  J 

The  functions  f  and  /.  will  be  assumed  to  have  con- 

tinuous  derivatives  with  respect  to  each  variable, 
a'  «  fQ(x',y')  is  the  position  of  the  "mean  camber  sur¬ 
face,"  and  2f  (x’t y'i  is  the  local  thickness.  Thus 
the  equation  of  the  hull  surface  can  be  written  in 
the  form 

(x'ty',M')  -  a'  -  [/  (x\y'i  t  f  (x'ty')]  *  0 
1  0  t  (1-1) 

When  the  body  is  moving,  its  trim,  heel  and  the 

position  of  its  center  of  gravity  relative  to  Oxyz 

will  change.  Let  a  be  the  trim  angle,  measured 

positively  in  the  bow-up  direction,  0  be  the  heel 

angle  measured  positively  clockwise  and  let  **  be  the 

amount  by  which  the  origin  O'  is  raised  (see  Figure 

2) .  Then  Oxyz  and  Oxyz  are  related  by  the  equations 

x'  »  *  ooea  +  [ (y  -  h)  oca®  +  a  ai«6]  si/to 
y'  m  -x  ainoi  +  l  (y  -  hi  cos®  +  z  ain®]  oo8<x 
a ’  m  -(y  -  hi  n«0  a  o oa0 


and 


x  *  x'  so 8 a  -  y '  aina 


(1-2) 


y  =  h  +  [x'  etna  +  y '  sob  a]  oob6  -  s'  sind 
b  *  [x'  eina  +  y'  ooea]  ein6  +  z'  oobQ 


Figure  2 
Coordinate  Axis 


Let  us  define 

H(x,y,z)  -  -(y  -  h)  einB  +  z  coe«  -  [f^x1  (x,yiZ),  yf(x,y,z)) 

+-  (I-3) 
which  describes  the  hull  in  the  system  Oxyz. 

The  motion  of  the  fluid,  which  is  assumed  irrota- 
tional,  is  most  easily  described  in  the  systems  dxyz 


/ 
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or  Oxyz .  Let  *(Z,y,B)  be  the  velocity  potential  in 
the  fixed  system  and  *(x3ytz)  that  in  the  moving 
system.  Then 

*(S3gtl)  «  *(x  -  ue  0086' tt  g3  S  -  u0  eini't)  (1-4) 

Both  the  velocity  potentials  4  and  $  satisfy 
Laplace's  equation 

5  *5*  +  *gg  +  *ziB  0 

5  ♦«  *  ^  *  ♦as  '  0  (I’5) 

The  absolute  velocity  of  the  water  described  in  the 

moving  coordinates  Oxyz  has  the  components 

(u  -  Uo  aoai ,  v,  w  -  u0  ain6)  «*  (6  -  u0  oob6,  *  .  4  -  u0  einfij 

x  y  a 

This  velocity  is  continuous  everywhere  in  the  fluid 

except  possibly  at  the  wake,  which  will  be  assumed 

to  occupy  a  thin  sheet  that  joins  the  body  at  its  aft 

end  and  extends  infinitely  along  the  negative  x-^xis. 

For  more  discussions  concerning  this  assumption  we 

refer  to  Thwaites  (1960)  or  Robinson  t  Laurmann  (1956, 

1.15).  The  pressure  in  the  fluid  can  be  computed 

from  Bernoulli '8  integral,  once  $  is  found,  as 

P/p  «  u0(6  0086  +  $  ain6)  -  gy  -  l/2(tl  +  +  **)  . 

*  (1-6) 

Let  the  equation  of  the  free  surface  be 
y  «  Y(xtz)  -  Y(x  -  u0oob6,  I  -  u08in6)  . 

(1-7) 

Then  the  potential  function  $(x3y,z)  must  satisfy  two 

boundary  conditions  on  this  surface.  The  first  is  the 

kinematic  boundary  condition 

*x{xt  Y(xtz)t  z)  Yx(x,z )  -  4  +  Yz  « 

*  uB(Y  oob6  +  Y  8in6)  . 

0  X  z 


(1-8) 
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The  second  is  the  dynamical  boundary  condition 

gY(x,z)  -  u0[$x(xt  Y(xta),  a )  gobs  +  ]  + 

+  +  by  +  ]  »  o  .  <I-9) 

There  are  corresponding  kinematic  and  dynamic 

boundary  conditions  to  be  satisfied  on  the  body’s 

wetted  surface  S  .  The  kinematic  condition  is 

W 

-(u  -  u0ooe6)  H  (x,y,z)  -  vH  (x,y,z)  -  (w  -  u0sin6)  H  (x,y,z)  «  0 
x  y  z 

♦Js  Hx(x,y,z)  +  ^|5  Hy(x,ytz)  +  ^|5  HJx,y,z)  = 

=  u0lH  (x,y,z)  ao8&  +  H  (x,y,z)  8in&]  ,  (1-10) 

where 

Hx(x,y,z)  *  -(fQx,  t  ftz,)  gob a  +  (f^,  i  f^,)  aim  , 

Hy(x,y,z)  *  -sine  -  [(f^,  t  f^,)  aim  +  (f^,  t  f  ,)  coaa )  coed 
Ha(x,y,z)  *  aoee  -  Uf QX,  t  ftx,)  aina  +  (f  ,  ±  ft  t)  gob  a]  sine 

Cy  (I_U) 

The  dynamical  boundary  conditions  on  the  body  are 
simply  the  equations  of  static  equilibrium  of  the 
forces  acting  on  the  body.  These  equations  can  be 
written  down  in  various  ways  to  conform  to  the  physical 
situation  of  the  problem.  Here  we  will  confine  our¬ 
selves  to  finding  the  hydrodynamic  force  components 
acting  upon  the  body  and  the  moments,  about  the  origin 
O'.  Denote  the  projection  of  the  wetted  hull  onto  the 

rm 

Oxyz  plane  by  swp  ,  the  water  pressure  on  the  starboard 
side  of  the  hull  by  Pg(xty,z)  and  the  pressure  on 
the  port  side  of  the  hull  by  P^(x,y,z)  .  Then 

?  -  [[  P(x,y,z)  n  da  , 
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1 

■> 

JL 

up 

0  0 

[p9(*,v,»)  fBX,  -pp(*>v,*)  fpx * )  * 

a 

5 

up 

(Parx,y,s;  Pp(xtyta)  fpv,]  dx'  dy'  t 

V- 

ils 

up 

[  -P8(x,y,a)  +  P^fx^,*))  dx'  dy’  . 

(1-12) 

s.fj 

P(xtyta)  'ft  X  n)  da  , 

0 

u 

vl 

t  i 

up 

pa  •  (-y'  -  -  yV  *  **'  dy •  , 

.  -f 

y  1 

'*< 

[ps  *  (z'  '  4*'  +  x>)  -  pp  ‘  r-a'^px'  ♦  dr'  ^ 

mx  •  1 

V 

[pa  *  V  ■  ^x''  “  PP  *  "V  •  *  V'}  * 

(1-13) 

A  kinematic  condition  muat  also  be  satisfied  on 
the  ocean  bottom.  Where  this  last  is  assumed  to  be 
of  infinite  depth,  the  condition  may  be  written  as 


lim  4 »  0  . 
y+~  ° 


(1-14) 


Finally  there  are  the  conditions  at  infinity, 
to  insure  that  waves  will  only  follow  the  body 

0([a:2  +  a2]  "l/2}  aa  xz  +  z2  •  for  x  >  0  , 

$(x,y,*)  «  ■  _  , 

0(1)  as  x2  +  az  *  for  x  <  0  . 

(1-15) 


Various  modifications  of  this  problem  are  possi¬ 
ble,  depending  upon  the  physical  situation.  In  parti¬ 
cular,  the  wave  resistance  of  a  catamaran  can  be 
determined  if  we  add  one  further  condition  analogous 
to  (1-14)  to  be  satisfied  on  its  plane  of  symmetry. 


II.  Method  of  Solution 


One  of  the  properties  of  the  problem  formulated 
in  the  preceding  chapter  which  makes  it  mathematically 
intractable  is  that  it  is  nonlinear.  In  order  to 
obtain  a  solution,  we  will  use  the  mothod  of  perturba¬ 
tion  expansion  as  an  approximation  to  linearize  the 
problem. 

Perturbation  Expansion: 

It  is  obvious  from  the  nature  of  the  problem 
under  consideration  that  the  disturbance  caused  near 
the  free  surface  is  dependent  on  two  perturbation 
parameters.  One  represents  the  "thinness"  effect  of 
the  body  8  ,  where  in  the  limit  as  8  0 ,  the  body 

degenerates  to  a  cambered  plane  of  zero  thickness. 

The  other  parameter  e  describes  the  asymmetrical 
effects  due  to  camber  and  incidence  angle.  As  e  -*•  C  , 
the  problem  becomes  one  of  a  symmetrical  body  aligned 
with  the  incident  flow. 

We  begin  by  imbedding  the  hull  form  (1-1)  in  a 
family  of  hulls  as  follows 

a'  *  t  3/W  (x’,y')  . 


Also  5  can  be  written  in  the  form 


(II-l) 
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6  »  z6  . 

(II-2) 

We  tentatively  assume  that  all  the  physical  variables 
can  now  be  expanded  as  an  asymptotic  series  in  terms 
of  the  two  parameters  6  and  z  .  Thus  we  obtain 
as  the  basic  expansions 

«  B$(lu)  (x,y,z)  +  (x.y,z)  +  Be*Ui> 

Y(x,z;8,z)  »  By^foa)  +  eyt01^  (x.z)  +  Bcy(n*  +..., 
a (&,z)  ■  Batl0)  +  Bea*11^  +..., 
h(&,z)  »  Bh<10)  +  Be/itn) 

era^e;  -  bb(10)  +  Bee(u'  +■••  • 


(IX—  3) 

Here  it  should  be  noted  that  this  analysis  may 
not  be  applicable  where  the  solution  of  the  problem 
is  singular ,  such  as  near  the  leading  edge  of  the  body. 
We  refer  to  Van  Dyke  (1964,  4.4)  for  more  discussion 
concerning  the  treatment  of  such  problems. 

The  problem  may  now  be  linearized  by  substitu¬ 
ting  the  foregoing  expansions  in  Laplace's  equation 
and  the  boundary  conditions,  besides  expanding  in  a 
Taylor  series  where  necessary,  and  collecting  terms  of 
the  same  order.  The  result  of  these  operations  is  a 
sequence  of  linear  boundary-value  problems  for  the 
potential  functions  4* ^ 1  ^ j  •••  •  The 
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first  of  these ,  which  represents  the  differential 
equations  for  the  potential  function  $(01)  ,  is  the 

well-known  problem  for  the  wave  resistance  of  thin 
ships.  The  velocity  potential  ,  which  is 

the  first-order  term  due  to  the  asymmetry  of  the  body, 
must  satisfy  the  following  equations  : 


A*<01>  =  0 


y  <  o  , 


*xzl)  (x>°>z)  *  K^(01)  -  0  >  K  "  9/“l  > 

Km)  (x,y,±0)  =  U0(6n)  -  fU)  )  on  S0  , 

"  ox 


lim  4> 5° (x,y,z)  *  0  , 


,101), 


0{[x2  +  a2]'1/2) 


0(1) 


ae  x2  +  z2  “  /or 


a:  >  0  , 
*  <  0  . 


(II-4) 


The  free  surface  y<01)  is  determined  by 


y(01)  (x,z)  =  Ko/^01)  (x,Otz)  . 

X 


( I I— 5 ) 

The  velocity  potential  <t>  ^  ^  represents  the  lowest- 
order  term  due  to  the  combined  effect  of  thickness 
and  camber  and/or  the  incidence  angle.  The  solution 
for  this  potential  function  will  not  be  given  here, 
though  it  may  be  of  value  for  investigations  in  the 
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future  .  We  will  confine  oureelvee  in  the  eubeequent 
analysis  to  finding  the  velocity  potential  ♦^#l 
First  let  us  apply  the  method  of  Green  functions  to 
solve  the  boundary-value  problem  (II-4)  for  this  velo¬ 
city  potential. 

Method  of  Green  Functions: 

As  is  usual  in  this  method,  it  relies  upon  the 
ability  to  construct  a  function  of  the  form 

G(x,y,a;  5,  r*  l+G0(xtyrBi  5,  n, 


where 

r-[(x-z)2+(y-n)2+{*-i)2)l/2 

and  G0  is  harmonic  in  the  region  occupied  by  fluid, 
or  in  the  case  of  this  linearized  problem,  in  the 
region  below  the  equilibrium  free  surface. 

Dimensionless  quantities  will  be  used  in  all  the 
expressions  hereafter,  by  comparing  lengths  to  half 
the  length  of  the  body  %  and  velocities  to  u0 
Also, 

l(xfyiB)-ueW($tyta),  G-l/lG  . 

Equations  (II-4)  will  then  become,  after  droppirg 
the  tildes, 

y<0> 

yo*8l/u\. 


*J°1}  (x»0,z)+yo$jj01)m0> 
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♦f01)  (x>V,±0)  -  «0)  -  on  S0  , 


Urn  $*ol)  (xty,z)  «  0  t 
y*-a>  y 

(o([r2  +  *2p/2; 


*  >  0  , 


<28  X2  +  z2  “  for 


[0(1) 


!x  <  0  . 


(II— 6) 

where  so  is  that  part  of  the  xy  plane  which  is  bounded 
by  -2  <  x  <  1  and  -T  <  y  *■  0  ,  and  T  in  the  dimen¬ 

sionless  draft. 

We  require  the  following  of  G  : 

-  o  , 

G^I'Xjj/jZyC^rU  +  Y0  Gn  *  0  , 


Ztm  G  «  0  , 
it*-«  n 


[0(U2  ♦  ;2]-1/2; 


ae  52  +  x2  -*■  ®  /0r 


5  <  0  , 


U  >  0  , 


(I  >7) 


Consider  now  the  region  of  fluid  bounded  by  the 


free-surface  plane  i  ,  the  two  sides  of  S0  and  the 
wake  w  ,  a  circular  cylinder  zR  with  Oy  as  the 
axis  and  a  radius  R  ,  and  a  horizontal  plane  at 

y  *  closing  the  bottom  of  the  cylinder.  Then  by 

Green's  Theorem  we  have  the  following  formula 
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4w4<01>(p;  -  |j  G(P;K* r\,0  -  4t01)^v)  *  * 

S(Q) 

-  \\  |(tc<Q1>  (t*,+0)  G(P;l*n>0)  - 

S0H) 

-  4<(U)  C?rP;C,n^;)  - 

-  (♦<°nfCJn,-W  crP;C,n,o;  - 

-  4*01>  (K,i\»"0)  ff?rP;C,n,o;)} 

♦  |J  (^01)  (i ,0,0  G(P;K,OtV  -  4(0l)Cr|)  dtfr\ 

lF 

-  Ztw 

n-K.< 

“B 

+  |2n  de  dn  P^01*  c  -  4 (01)  gr)  . 


(^0l)  0(PsltntO  -  ♦{01^n) 


Substituting  (II-6)  and  (II-7)  in  the  above  formula, 
we  obtain 


4jt+M  (P)  -  -  ||  G^fPjij^OJ  [4><01>  fC,n,+^  -  $l01)  dedn 

-  j|  -  *(01^  «fc«fc 

f2*  [0 

■#•  d0  dn  0(1)  . 

J0 

The  seoond  into ;:r:..l  :nry  be  integrated  by  parts  to  give 


I 


where  this  line  integral  is  taken  around  the  interaction 
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of  Zp  and  .  As  R*»  ,  both  this  integral 

and  the  third  one  converge  to  sero.  We  finally  have 


*<01)  (x. 


G^(P;Z,t),0). 


{Ztn*-0J)d£dr\ 

(II-8) 

It  should  be  noted  here  that  the  above  analysis 
was  done  in  more  detail  following  Wehausen  (1963)  and 
the  same  final  result  was  obtained. 

The  solution  for  the  Green  functions  is  well 
known  for  several  physical  situations,  many  of  which 
are  given  in  Wehausen  and  Laitone  (1960) .  Here  we 
use  one  form  given  by  Eggers  and  others  (1960) ,  which 
is  suitable  in  our  analysis. 


[  K((y+r\)+iw  ) 


K-yeec^Q 

O 


dk+ 


1  dk 


Kc-y  sec20 


(  TT—Q \ 


where  and  a ^  represent  the  indented  paths  of  the 
contour  integration  around  the  pole  as  shown  in  the 
figure  below 
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i 


| 


1 

\ 


and 

r^((x-i)z+(y+r\)2+(z-t,)2)  1/2i 

i  M(x-Z)co8Q+(z-t;)8ine, 
l 

In  handling  the  above  countour  integrals,  we  should 
be  cautious  since  its  integrands  are  not  only  oscilla¬ 
tory,  but  also  possess  a  singularity  at  k-y»  aeo2e 
for  each  value  of  e  .  It  turns  out,  however  (see 
Appendix  A) ,  that  by  a  proper  choice  of  integration 
path  in  the  complex  plane,  this  integral  can  be  simpli¬ 
fied  considerably,  namely 


o 
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where 


p  *y  a#<?20  l(y+r\)+iQ  1,  p  *ye  8eczQ[(y+r\)-iiL  ]  . 

1  12  1 

and  H(x)  is  the  Heavyside  function,  defined  by  H(x)~l 
for  x  >  0  and  H(x)*0  for  x  <  0 

Differentiating  this  function  with  respect  to 
C  ,  we  obtain 


G^(x,y,z;^,r\,0) 


Re  lo 
it 


ae<?40  ain0  d0- 


X  i\ 


e  P[ - I - 1 - 1  dp  - 

o  [('p+p^2  <'p+p2-’2J 

f°  y  *2(y+n) 

e  0  ein(y 0\(x-V ). 


X  ain[y0\  -  1  *  a)  A*  d\ 
f°° 

-2y  6(-ui)  einl2y  \(x-£)]  d\,  (II  -  11) 

O  1  O 

where  we  have  made  the  substitution  \  *  eec6  in  the 
last  integrals. 


Determination  of  the  Volocity  Potential  a*01) 

The  integral  equation  for  the  velocity  potential 
,  as  in  (II-8) ,  represents  a  distribution  of 
doublets  over  the  surface  sg  and  the  wake  w  .  The 
moment  \i(x,y)  of  such  a  distribution  is 

4i\i(x,y)  •  (x,y,+0)  -  (x,y,-0). 

(11-12) 

Therefore,  the  velocity  potential  $(ui)  ,  after  deleting 
the  superscript  <01)  ,  can  be  written  in  the  following 


form 


(11-13) 


4 (x,y,a)  «  Jj  -G^(x,yta;S.tT\tO)  didr\  . 

S0+u 

Now  the  linearized  dimensionless  pressure  (based 

on  the  pressure  i.  p «02  )  associated  with  this  velo- 

2 

city  potential  can  be  found  from  (1-6)  as 
P(x,y,z)  «  2$x(xty,z) 

Accordingly,  the  pressure  jump  across  the  plane  a  *  0 
is 

P(x,y,+0)  -  P(x,y,-0)  -  2  Ux(x,y,+0)  -  *x(x,y,-0)]  , 

«  2  l$(x,y,+0)  -  $(x,yt-0)]x  . 

From  (11-12) ,  this  becomes 

P(x,yt+0)  -  P(x,y,-0)  -  8V]xx(xty) 

(11-14) 

Since  the  pressure  must  be  continuous  across  the 

wake,  then  from  (11-14),  u  must  vanish  on  the 

x 

wake.  It  follows  that  y  is  constant  along  lines 
parallel  to  the  x-axis  in  the  wake  and  therefore 

\i(i }y)  =  v(T,y)  ,  for  x  <  T  . 

(11-15) 

where  *  *  T(y)  is  the  equation  of  the  aft  end  of  the 
body.  Also,  since  the  velocity  potential  must  be  con¬ 
tinuous  everywhere  outside  the  surface  s0  and  w  , 
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therefore  it  follows  that 


\i(x,-t)  =  0 


(11-16) 


On  the  other  hand,  due  to  the  expected  singularity  at 
the  leading  edge  of  the  body,  there  will  be  a  discon¬ 
tinuity  in  yx  at  this  edge  (see  Appendix  D) . 

Equation  (11-13)  may  now  be  written  as 


$(xty,z) 


-G  (x,y,z;t,r\t0)  v(Tty)  df.dn  + 
"w 


+ 


f  G  (x,y,z;t,,c\,C) 


didr\  . 


(11-17) 


The  velocity  potential  $  can  be  then  determined 
once  the  distribution  function  v(x,y)  is  known. 

To  find  y  ,  we  apply  the  linearized,  kinematical 
boundary  condition  ( II— 6 )  on  the  surface  S0  to  the 
above  equation: 


lim  3 
z+0  33 


G  d^dr\  + 

.  W 


G  \i(l,r\)  d^dr\ 


-  +  fjx>y} 

ox 


(11-18) 


This  is  a  Fredholm  integral  equation  of  the  first  kind 
for  the  unknown  doublet  moment  u(x,y)  .  The  closed 
analytical  solution  of  this  integral  equation  is  beyond 
our  resources  because  of  the  complicated  nature  of  the 
kernel  function.  Consequently,  a  numerical  approach, 
based  on  the  finite-element  method,  will  be  used  in 
the  following  chapter  to  obtain  an  ■  pproximate  solution 
for  the  doublet  moment  v(xty)  . 
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III.  Solution  of  the  Integral  Equation 

For  convenience  in  subsequent  calculations,  the 
surface  S0  win  be  approximated  by  a  rectangular  plane 
and  the  coordinate  system  Oxyz  will  be  shifted  to  the 
system  shown  in  Figure  3.  Thus  T(y)  *  0  .  Let 

us  divide  S0  into  rectangular  elements  determined 
by  a  finite  number  of  nodal  points  as  in  Figure  3. 


y 


The  Surface  S0  Divided  into  Finite  Elements 


Now  in  order  to  solve  the  integral  equation  (11-18) 
approximately,  but  to  any  required  degree  of  accuracy, 
it  is  convenient  to  assume  that  the  presumed  piece-wise 
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continuous  function  \i(xty)  over  the  surface  S0  can 
be  approximately  determined  by  a  finite  number  of 
values  at  the  nodal  points,  as  described 

below.  Where  the  nodal  values  of  u  define  an 
approximating  function  within  each  element,  the  integrals 
in  (11-18)  can  be  performed  analytically  over  the 
elements  and  over  the  semi-infinite  strips. 


y 


Figure  4 

A  Rectangular  Element 

Consider  the  rectangular  element  shown  in  Figure  4. 
If  the  function  u  is  to  be  continuous  between 
adjacent  elements  it  is  necessary  for  u  to  vary  in 
a  linear  way  along  the  sides  of  the  elements.  Then, 
coincidence  of  w  at  the  nodes  will  automatically 
Insure  coincidence  at  intermediate  points.  The  function 
u  within  the  element  may  be  written  as 


-2*- 


\i(i,r\)  -  1 

4LiLr\ 


(Li  -  i)(Lr\  -  V  + 

ffljK 


+  (Li  +  i)(Lr\  -  nJ  +  (Li  -  i)(Ln  +  n )  v  .  + 

m+itn  m,n+2 


+  (Li  +  i)(Lr\  +  n;  V, 


m+ltn+l 


or 

v^n;  »W2,n  - 

4LiLr\  v. 

-  rWc;Vn;  y«+J 

for 


«m  <  «  ‘  Vl  '  \  ‘  Vl  ' 

(III-l) 

Now  by  satisfying  equation  (11-18)  at  a  number 
of  points  on  Sg  equal  to  the  unknown  p's  ,  the 
integral  equation  reduces  to  the  following  set  of 
linear  algebraic  equations  in  u'e  : 


N 

nil 


lim  3 
z-*Q  Tz 


fn”^  \i(0, nJ  dr\ 


G  (Xityj,z;itT\t°)  di  + 


N  M 

+  I  I  lim  L, 

n*l  ml  zM)  3* 


n”+2  dr\  G.(x.,y.,z;i,r\,0)  \i(itn)  didr) 

’  t  J 

5. 


m 


-  -«  +  4/*;^  • 

(III-2) 

To  facilitate  evaluation  of  the  above  integrals, 


we  will  define  the  following  variables: 
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t  -  x  -  e  , 


q  -  y  -  n 


a  ■  y  +  n  -  2t  , 


x(p,*)  *  p  +  (y  +  n  -  2x)  . 

'  A 


Then  G  may  be  written  in  the  form 


( III— 3 ) 


G^(t,q,ata)  «  »/(  t2  +  q2  +  a2)*/2  -  a/(  t2  +  a2  +  a 2)^2  -  Ra 


X  -  einQ  d6  e~p  [l/(  x  +  tw;2  -  l/(  x  -  iW2]  dp 

.  “t  Jo 


ffr  uj  eY 

*  1 


TTi  eY°X  8  ainTYoi  * 


*  ain(y *  aj  X2  dX 


-2y0  I"  6(  uj  «YoX*a  ain(ty0\t)  d\,  (III-4) 


where 


b>  «  t  ooaB  +  a  aine  , 


ii)  m  -  ^.(t  +  •  a)  . 


Also  equation  (III-l)  may  be  represented  as 

u(t,  (*  J)  «i‘I  ul  -  «?  \)  , 


(III-5) 


where  in  this  system  of  notations  (  ®  )  means  the  same 
equation  is  applicable  to  both  the  variables,  e  corres¬ 
ponds  to  the  upper  subscript  and  q  to  the  lower  one; 

A  is  the  r rer  of  the  element  which  is  equal  to 


,  [Hj  is  the  row  matrix 
[y]  *  [Vmtn*  Vm+ltn*  vmtn+l*  vm+ltn+l] 


( III— 6) 


=*!»*. 
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and  W  j)  is  ths  column  matrix 


{Q  i>- 


m+1  ~  t)  '  ((  8  )n4.i  ~ 


n+1 


+-(*m  -  *  "  q  >n+l  ~  (q  )} 


lit.  -  t) 
m 


Oq'n-'q" 


• 

( 

Vn+1 

ft m+1 

±( 

8 

<7 

}n+l 

« 

~*m 

( 

8 

<7 

}n+l 

It 

m 

+  ( 

8 

<7 

*n+l 

1 1 

ft m+1 

• 

( 

8  ) 

q  n 

+  ( 

8 

<7 

>n 

tl 

K' 

( 

8 

<7 

>n 

ifm 

t( 

8 

<7 

>n 

11 

1 

(  8  ) 

1  <7 

t 

,  • 

It  follows  that  (III— 5)  may  also  be  expressed  in  the  form 


i 

<°> 


where  [  (  j  )  ]  is  the  row  matrix 

<l>2.  <V,‘  'SJJ> 

This  may  also  be  written  as 


(III— 7) 


a  J]  *  i  '  [  v] 
B  A 


'Hml 

•  ( 

“ml 

±( 

tt  • 
m 

'? 

^  n+1 

+t 

m 

+( 

“ml 

’  ( 

8 ; 
q  n 

+* m+1 

*( 

U  * 

.  m 

(  8 
<7 

^n 

t( 

q  'n+1 


q  'n 


(III— 8) 

Consider  now  the  contribution  of  the  wake  to  the 
integral  in  equation  (III-2) ,  which  we  will  denote 


-25- 


ty  Jw(xi*V^  •  Hire 

»(*.  <’<,  J)  ■  <  S\  +  <‘q  > 


(III— 9) 


where 


r  0  )12  "  ZSn  ’  vl,n+l)  * 


*  *J 


(III-10) 


V-’e  may  then  write  Jw‘  as  follows: 

Jun(xi>y, j> " liM  fj  If!"2 

v  9+n  uHn 


z+0  W» 


AT  (a  dt/(  t2  +  q2  +  z2)3^2  )  +  Jfl  (o^  •*  a^;  de 
X  "  dt  { z/(  t2  +  a2  +  z 2)3^2  +  Re  \  ai«e  de 
X  |  e  p  tJ/r  X  +  iw;2  -  i/(  x  -  iu)2  ]  do  +  4y*0 


*oo  ^  2 

x  Jj  eYo  8  (y^einrY^^einry^A*  -  T  *  a;  *  A2  * 

+  (l/2yg)6(G)8in(2y\t)]  d\  .  (m-li) 

This  function  is  evaluated  in  Appendix  3,  and  in  re¬ 
ference  to  this  may  be  written  as 

bl,l  bl,2  bl,3 

J»n(xi>yj}  *  [8  U  62|  ’  '  '  +  [aj,  a2]  •  • 

■b2jlJ  b2t2  b2,2. 


This  may  now  be  expressed  in  the  form 


*V  (xi*y3}  "  1  VI, n»  vitn+l  J ' 


(jn,n.  .1 
\Bi  (xi>y.i> 


W2*n(xoVj)i 


( II I— 12 ) 


I 
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where,  from  (III-10) 


1)  b, 


\8n+l  -1\\bl,i  +  b1.3 


\»\,n  (xi w 


qn 


'“an  »  ^b2t2  *  h2,3> 
(III-13) 


Consider  next  the  contribution  of  the  surface 
to  the  integral  in  equation  (III-2) ,  which  we  will 
denote  by  Jn»n  (x.t\ j>)  .  From  (III-4)  and  (III-7) 

8  t  J 

this  may  be  written  as  follows; 

f  K+l  *ml  (&i  +  „ 

J”l,n  -  Um  i-i*  .  .  ,  ,  dt 


a  1  ,o0 


z+0  3s  [  }qn 


(t2  +  qz  +  z2)^2 


Bn+1  ,  Wj  r/  ,  ,  ... 

+  da  [(<*!  +  a2a;-f('o3  +  a4e;t]* 


(t2  +  a2  +  z2)^ 


■f  1  1^  ainQ  dfi  Re  i  e  [" - - - -  -  - - - rj  dP  + 

IT  J.|  Jo  |/x  +  iv>r  (X  -  tu>;zJ 

+  4'i2^^°^lB[ji(^)  ein(y0H)ain(y  xA 2-l‘z)'\2  + 

+  (l/2y0)6(G)8in(2y0\t))  dxj  dt^j  I 
It  turns  out  (see  Appendix  C)  that  this  expression 


(III— 14 ) 


may  be  written  as 

(x^yj)  =  [*K)'{aKtl)  +  l«J  •<«„,*  *  +  a<,4+aK,S} 

and  k  *  1,2, 3,4 

Substituting  for  (  ®  from  (III-8) ,  we  may 

express  d^*n  (x^,y^)  in  the  form 


(xityj)  *  ly] 


'  AT  <H>yj} 
Afn  (xi>yj} 

‘  Afn  (x.,y.) 

A>n  (*i>yj)' 


(III-15) 
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whara 


Afn 

(x.ty.) 

* 

t m+l 

qn+l 

~*m+l 

“Qn+l 

1 

*  ' 

al,l 

An,n 

A2 

1  . 

Qn+l 

Qn+l 

-1 

• 

a2tl 

Am„>n 

A3 

<xi-Vj> 

*  A 

mti m+l 

qn 

tm+l 

qn 

-1 

« 

°3,1 

Am,n 

A4 

<xi‘«j>- 

t 

m 

qn 

~lm 

-qn 

l 

■a4. 1- 

+1 

A 


*m+i  8n+l  i,n+l  8*+l 

al,2  *  al, 3  *  al,4  *  altS 

*m  8n+l  -tm  ~8n+l  1 

a2,2  +  a2,3  +  a2,4  +a2tS 

t  .  b  -t  „  -a  7 

m+l  n  m+l  n 

• 

a3t2  +  a3,3  +  a3t4  +  a3fS 

-t  8  t  B  -1 

m  •  n  m  n  J 

a4,2  +  a4,3  +  a4,4  +  <*4  c. 

(III-16) 

We  may  now  write  equation  (III-2)  in  the  form 


?  f4"  Mj'l 

»  ft 


5  2  (xi W 


m«2 


4  -  4*  te<*V  ■ 


ic  -  lt2,3t4 


(III-17) 

The  left-hand  side  of  this  equation  can  be  written  as 


{Sl’”  (xi-Vj>  H ,n  *  4’"  * 


*  J2  y. +  aT  'W  \n+i) + 

M+l  ,  * 

*  k=2  ^  *”  y<,«  4  (xi*yj}  yK,nv-l)} 

{[  V  +  4i,n  (H>y^  1  *  4?”  (xi*yj}  Vm,n  + 
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*  ?  w  (x^i  *  t?1-”  vy  v4 


+  I  { H'*'1  fvV  *  'vV  v +  ',"'K'!  ».«* 

K*o  V 

*  ?,  *  ATht'1  %.«} 

m*2  i 

!,  (  Bi,n  tefV  *  *j*"  'VV  “1,»  +  (xv»j>  *m j.»}  + 

*  2  Bl’n'1(x-i-vi>  +  isn'1,zi‘!lJ>  ui»n*  ^'"'S ^m.n}  * 

*  L  {  ',3,i  teiV  +  ,l2’1'3  HI  * 


+  <”■*  !xi>y^  *  wi^K+iJ + 

+  f  !.  + <r3-n  + 


*  /q-"-1  (xilVjl  *  %1>n-1  (x^y.)  „m<„ 

From  equation  (II-6)  ,  we  find  that 


um,  1 


0  , 


m  =  1,2,... ,M+1  . 


Finally,  we  can  write  the  set  of  the  linear  al¬ 
gebraic  equations  for  um  n  in  the  form 

N+l  M+l 

nl2  nil  Cm>*(Xi,y0)  V*  =  fox(xi»»j)  "  6 


(III-18) 


where 


m  i?nixi‘>jl  *  At1,n,xi-»3>  * 


.-29 


^r1  w  ■ 


m  »  2,3,...  ,M  ,  n  a  2,3,. . .  ,N  ; 


m  »  2,2,... ,M  ; 


cM+i,n(xi»yj}  *  A2,n(xiyy^  *  AM4n’1(xi>yj*  > 


n  »  2,3,...  ,N  ; 


VVV  "  Bl‘n,xi‘yj>  *  B2n'1<xi-yj>  * 


*  *  *Ii‘n'1<xi.yj>  , 


n  «  2,3,...  ,N  ; 
cM+i,N+i(xi>yj)  mAl4N(xi>Vj}  ; 

ci,m(xi*yj}  ■  Bl’N(xi>yj}  +  Al,llfxi*Vj)  • 

(III-19) 

Equation  (III-18)  can  now  be  solved  numerically 
for  a  given  body  moving  with  a  constant  velocity  ug 
at  an  angle  of  attack  5  . 
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IV.  Numerical  Results  and  Conclusions 


A  computer  program  was  developed,  based  on  the 
numerical  analysis  given  in  Appendix  D,  to  compute 
the  moment  distribution  y  as  well  as  the  forces  and 
moments  acting  on  a  yawed  and/or  cambered  body. 

The  numerical  scheme  was  first  applied  to  the 
case  of  zero  speed,  where  the  free  surface  is  regarded 
as  a  rigid  wall  and  th  s  flow  is  the  same  as  that 
around  a  fully  submerged  double  body  consisting  of 
the  ship  and  its  image  over  the  free  surface.  The 
curves  in  Figures  5  and  6  show  the  results  of  such 
computations  and  are  compared  to  similar  results  given 
by  Thwaites  (1960,  p.  343)  for  rectangular  wings. 

The  effect  of  the  free  surface  was  then  considered 
for  a  small  Froude  number  (Fn=  u0//gL  =  0.1)  and  the 
computations  were  made  for  two  different  values  of 
the  length/draft  ratio,  L/T- 7  and  L/T=20.  The  results 
showed  a  very  slight  deviation  from  the  values  obtained 
for  the  case  of  a  rigid  surface. 

Computations  were  carried  out  for  a  yawed  body 
having  a  length/draft  ratio  equal  to  20  at  different 
Froude  numbers  and  two  values  were  used  for  the  number 
of  nodal  points  in  the  longitudinal  direction,  M- 5  and 
M<*  10.  Also  one  run  was  made  for  the  same  body  with  M*  20 
and  Fn»0.34.  The  results  of  these  computations  are 
shown  in  Figures  7-11. 
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There  is  a  large  discrepancy  between  the 
results  obtained  when  M  is  set  equal  to  10  and  those 
obtained  for  W* 5  and  A/*20  in  the  range  of  Froude 
numbers  between  Fn»0.225  and  Fn«0.45.  Norrbin  (1960, 
p.  379)  pointed  out  that  "for  a  surface  ship  running 
at  Froude  numbers  exceeding  Fn  =  u0//gL  -  0.3,  wave 
formation  is  generally  found  to  cause  a  change  in 
trim  and  stability  characteristics."  Hu  (1961) 
found  out  that  the  lateral  stability  derivatives  have 
increased  about  50%  above  their  values  at  zero  speed 
for  a  Froude  number  Fn»0.35.  Due  to  the  lack  of 
experimental  daca  at  higher  Froude  numbers,  it  is 
difficult  to  predict  the  behavior  of  the  force  and 
moment  coefficients  in  this  region.  On  the  other 
hand,  there  is  the  possibility  that  the  numerical 
analysis  is  divergent  for  the  particular  value  Af-10 
in  the  range  of  Froude  numbers  where  this  behavior  is 
encountered. 


Conclusions 

Although  one  of  the  purposes  of  this  study  has 
been  to  investigate  theoretically  the  side  forces  and 
moments  acting  on  a  yawed  body  in  a  free  surface,  an 
equally  important  one  has  been  the  testing  of  the 
applicability  of  a  numerical  method  for  solving  the 
complicated  integral  equation  that  arises.  The  agreement 
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of  the  results  obtained  by  this  method  with  those 
obtained  by  other  methods  when  the  Froude  number  is 
small  seems  to  indicate  that  the  method  is 
fundamentally  sound.  For  higher  Froude  numbers  the 
inconsistency  of  the  values  tor  Nm  10  with  those  for 
H»S  and  W-20  indicates  that  further  investigation  of 
numerical  stability  in  this  region  is  necessary.  Further 
experimental  work  supplementing  Norrbin's  measurements 
is  also  desirable. 
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LIFT  COEFFICIENT  ON  RECTANGULAR  WINGS 
No  Free  Surface 
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Figure  5 


Figure  6 
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Figure  7 
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Figure  8 
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ADDED  RESISTANCE  COEFFICIENT 
For  6  deg.  Yaw  Angle  L/T  *  20 


Fn 


Figure  9 


Cmx/6 


HEELING  MOMENT  COEFFICIENT 
For  Yawed  Body  L/T  *  20 


Figure  10 
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Figure  11 
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Appendix  A 


Green  function  G(x,ytz;Zt is  given  by  equa¬ 
tion  (IX— 9 ) .  Consider  now  the  integral  with  respect 
to  k  in  this  function,  which  we  will  write  in  the  form 


d< 


<  -  y oSea^Q 


j-  exp  [-<ae  ^  ]  + 


ic  -  Y  8ec*& 
° 2  0 


d*  exp  [-< ae^  ] 


(A-l) 


where 


-ae  *  (y  +  r\)  +  iu  , 


a  >  0  . 


Therefore  -a  aoep  =  y  +  n  <  0  , 

2  ~  ~  2 

and  u  =  a  sinp  . 

Let  us  evaluate  the  above  integrals  over  the 
straight  line  < =ret^°  in  the  complex  k  -  plane.  Then 
we  obtain 


- exp  \rarei(^0'^  Jf 


0  rel*#-Y  0sea2Q 


dreUo 

re^°-Y  .sec2  9 


exp  [~ave^^0+^>1  ] 


The  oscilatting  part  in  these  integrals  will  vanish 
if  we  set  p  =  p  in  the  first  one  and  P0=  -<l»  in  the 
second.  Then  it  becomes 


,-ar 


jt  r  -  y 08ea*Qe 


a>  r  -  y08ec*Ce 
2 


, far  .  dr  = 

- D* 


,-P 


o'  P  *  P. 
"i  1 


dp  + 


9-P 


p  +  p 
°2  2 


dp 


(A-2) 


pj  *  y g8eo26  [(y  +  r\)  +  iw] 
P2  *  y0seo2Q  i(y  +  r\)  -  tu>] 


where 
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We  will  apply  now  Cauchy's  Theorem  for  contour 
integration  to  evaluate  the  integrals  in  (A-l)  .  The 
paths  of  integration  will  be  completed  as  shown  in 
the  figures  below  to  form  a  closed  contour. 


For 


Gi 


<  o 


In  order  that  the  integrals  along  the  arc  converge 
to  zero  as  R  -*■  •  ,  it  is  necessary  and  sufficient  that 

-tt/2  <_  4>0  t  iji  <_  w/2  ,  This  requirement  is  satisfied 

by  the  previous  choice  of  $0  . 

If  <3  <  0  ,  we  must  take  into  account  the  residue 

at  the  pole  at  <  ■  y oaec2e;  if  2  >  0  ,  there  is  no 

residue.  By  using  the  result  given  in  (A-2) ,  the  integral 
in  (A-l)  can  be  written  as 

"  e"P  f,  1  *  1  .  1  dp  +  H(-G)  {Residue) 

[mr  mrj 

JO 

where  H(-&)  is  the  Heavy  side  function. 
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Finally  we  will  now  evaluate  the  residue  at 
<  *  Y<>  sec2  0  as  follows: 


-Re  IS. 

IT 


$  dQ  sec2*  H  [-e'o  89a'*  [  <«  +  *)  +  i»  ]  + 
+  ey0  secz6  [(y  +  t))  -  ioj]  I 


Im  y0 


.a. 

2 


fT 


Im  2ye 


-  -4yc 


de  ecch  [ey°  Se°2e  1  *  ”>*  - 

^y o  seo2e  [(y  +  r\)  -  tui]J  + 

i,  We  [>”  SesH  ltt  *  ti“1 

y0  eea!e  \(y  +  nJ  -  iu]j 

x  , 

de  seo2d  eYo  eea  0  ^  +  n>>  *  coe[y0(z  -  r.)  ain%  aec20] 

0 

X  [giy O  sees  (x  -  V  g-iy 0  seed  (x  -  £)] 

x 

ey0  8ea  B  (y  +  r\)  ,  g^n[y^  egflg  (x  _  ] 

0 

X  ooe[y0(z  -  Z)  sinB  seo2B  ]  sec20  d6  (A-3) 
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The 


Appendix  B 
jfj*  Integral 


This  integral  is  defined  by  equation  (III-ll) . 
We  may  write  it  in  the  following  form; 

Ii  -  f"+3  <\  *  \«>  * 


X  P  dt/(t2  +  q2  +  z2)^2 


' xi 


-z 


n+1 


fa  +  a  e)  da  |  dt/(t 2  +  a2  +  z2)  ^2 

1  2  lx. 


,xi 


-  ~  |  dt  \n+1  fa^  +  ct^s)  ds  j?  sine 

X •  0  ' 


de 


t  n 


r 


*  1  a'°  '  s*  <7x  +  TZ!*  -  TxTTZF1  *  * 


*V,  J”  it  j  n+J  fa  *  *  |”  «Y“X‘*( 


’X.  '8 

i  n 


[2yQH(  a) ; 


AT  ainfYoAt;«£n(YoX/\s-^;*X2  +  6(  u)ein(2y  \t))  d\ 


(B-l) 


Now  we  will  evaluate  each  of  the  above  integrals 
in  (B-l)  separately  and  determine  the  limit  of  its 
derivative  with  respect  to  z  as  z  +  0  .  Consider 
the  first  integral,  which  we  will  write  as 

-  -  llT  »,  *  V  *  t  7777777* 
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After  integrating  this  once  over  t,  we  have 


The  integration  of  the  second  term  in  each  of  the  above 
integrals  can  be  easily  performed  if  we  make  the  sub¬ 
stitutions  Vj  * 
in  the  second.  Thus 


The  contribution  of  ^  to  is 


in  the  first  and 


i 


(B-3) 


Mmavnur-*. 


where 


■  LI  L 

>1‘1  ’  \%+l  '  In)  '  H  [Vl TTj* 

%i  *  Zn(raJ '  "**rl^  ■ 


ooth~l\-£-^j  -  ooth~l  ~yt 


(B-4) 


The  second  integral  in  (3-1)  may  be  written  as 
[8n+l  r 


J2  «  -3 


(a  +  a  e)  de 
1  2 


dt 

(t2  +  8*  +  Z^j- 


By  compariso”.  /ith  Jj  this  rr  y  now  be  expressed  as 

r  i  8n  j  8n+l  ] 

h  -  *x  ■  r- tan  ~r 

*  [of1  M  Txfl 

+  *  32J  -  ln(8n+j  +  32jJ  - 

-  3  [oath''  **  *  *  +  ""+1  -  ooth-'  /Z^X".:,an 

L  H  H 

The  contribution  of  this  to  ^  (xi*Vj)  is 


2  >  + 


1'  r  i  . 

Ztm  _LL  *  ( a  ,  a  ] 

Z+0  12  br,  n 


(B-5) 


where 


'  1 

1  ■ 

i 

sUi 

4J 

.5n77  " 

K 

'  H 

8 n+1 

en  . 

In  JjL 


ooth~ 1 


st,J 

-  1  wr~ 


(B-6) 
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In  the  third  term  in  (B-l) ,  the  integral  over  t 
has  an  improper  upper  limit  and  consequently  changing 
the  order  of  integration  would  not  be  desirable.  To 
avoid  this  difficulty,  we  write  this  term  as  follows, 


t*/2 

)de  sin0  dB  . 

-n/ 2 


R 


e 


i 

( X  +  iu) 1 


i _ 

?  X  -  ■£<*>;* 


]  dp  . 


Now  this  integral  is  in  a  proper  form  and  is  similar 
to  the  integral  I3  in  Appendix  C.  Therefore  the 
contribution  of  this  to  J™*n (x^,y j)  can  be  expressed 
in  the  form 


lim  ^3 
a+0  3s 


I  )  • 


(B-7) 


where 


b 


1,3 


9 


b 


2,3 


k 

Z 

k-i 


a 


2,3 


(B-8) 


For  numerical  computations,  the  upper  limit  of  summation 
has  changed  to  a  finite  value  k  that  is  found  to  achieve 
the  desired  numerical  accuracy. 


Finally  we  will  evaluate  the  contribution  of 
the  last  term  in  (B-l) .  The  value  of  this  integral 
will  depend  upon  the  sign  of  u>.  Since 


0  «  (x.  -  cos0  +  z  sine 

Z* 


where 


K  <  0  and  x  .  >  0  , 


therefore  the  sign  of  13  will  depend  upon  z  sin0. 

In  the  limit  as  z  -*■  0,  35  >  0  and  consequently  H(-iL)"  6(-Ui) 

Hence  it  follows  that  the  contribution  of  this  integral 
j!  yt 

to  Ju*  (xityj)  vanishes.  This  result  is  expected  since 
this  integral  represents  the  contribution  due  to  the 
waves  generated  behind  the  doublets  at  £  »  0.  We  now 


have  the  following  formula  for  j2*n(x,,y„ 

w  t  j 


1 

*  I°7>  ' 

2,1  2  2 


b l,2*b  1,  3 
b2,2+b2,3d 


(B-9) 
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APPENDIX  C 


The  (xj*V Integral 


This  Integral  is  defined  by  equation  (III-14).  Applying  the  formula 
for  integration  by  parts  to  the  third  term  of  this  integral,  we  obtain 


_2£ 


yQ  aeo  0  Re  i 


III)  x  -  10) 


"  Y0  80O2B  |  -  J 

\  8  +  (t)  J 


-p  2o> 

*22 
x  +  w 


dp 


-  yQ  aeo°B 


15 

j 


,-P 


2d) 


2u) 


2  2  2  2, 

8  +  U)  X  f  111  J 


dp  . 


After  changing  the  order  of  integration  of  the  last  three  integrals  in 
®*y  be  written  as 


-  Um  | 


8  t 
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5n+l 

J  dq 
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TT 

Y  F  ^n+l  ~ml 
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J  IT  J0  J  a  J  /■»  l  1  2 


a  w 
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J>n+1  2  r 
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J t  H(u)  •  ein( yo  XfcJ  • 
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12  SO 


+  2yo  j"  <*A  j*n*J  «YoX2®  da  jl”*1  6(  u)  ain(2y0w  • 

•  [fct^  *  +  (a^  *  a  a)  •  tj  dt  (c-l) 

Each  of  tha  above  integrali  will  be  evaluated  separately;  also  the  limit 
of  its  derivative  with  respect  to  a  as  a  *  o  will  be  calculated.  Consider 
the  first  Integral,  which  can  be  written  as 


rn+1  fml  (&.+  S„  q)  +  re,  +  6 A  q)  t 

xi  m  *  L  *  Jt  ~:'.f  .  3*.J  * . 


(t2  +  q2  +  a2)3/2 
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Integrating  one*  over  t. 


A  similar  Integral  to  the  first  one  above  has  been  performed  In  Appendix  B; 
therefore 
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Differentiating  this  aquation  with  raspect  to  a  and  taking  the  limit  as  «  0 

wa  obtain 
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where 
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fora 


-57- 


lim 

*+o 


&2*  hJ 


(  \ 

au 

a3tl 

v  / 


(C-A) 


Consider  the  second  integral  in  (C-l),  This  can  be  written  as 
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Consider  the  third  integral  in  (C-l) ,  which  will  be  written  ae 
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By  integrating  the  laet  term  over  u,  we  have 
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Since  this  equation  is  In  a  proper  £orm,  then  It  le  possible  now  to  take 
the  limit  o£  the  derivative  with  respect  to  Z  as  Z-*o  before  integrating 
over  5.  Than  it  becomes 
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Integrating  now  over  8 1  we  obtain 
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(C-9) 


[a3,  o2>  a3,  a4] 


(C-10) 


-62- 


whare  froa  (C-6)  and  (C-7) 
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wh«r«  we  have  made  the  eubetitution  X  -  eeo9  in  the  above  integrale. 
Consider  next  the  fourth  term  in  (  C-X) 


Bn+1  „  ,2 


Bin  (y\  / X2’2  •  z;*X2<JX  J  e  0  da  f 
*  ~ 


H(u)  ain(y\t)  . 
o 


n 


n 


[(xx2  +  a 2  $)  +  (a3  *  a4  «;•  t]  it 


-64- 


It  is  possible  here  to  take  the  limit  of  the  derivative  with  respect  to 
*  as  s  ■+  5  before  integrating  over  t  and  e.  Then 
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Integrating  over  8  ,  we  obtain 
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Th«  contribution  of  this  integral  to  can  be 
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(C-15) 


Finally  the  last  integral  in  (C-l)  can  be  written  as 
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integrating  over  t  we  obtain 
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Differentiating  with  regard  to  a  and  taking  the  limit  as  s  -*■  0, 
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we  obtain 
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The  contribution  of  thia  integral  to  J^*n(x^ty^)  can  ba 
written  as 
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Appendix  D 

Numerical  Analysis: 

To  determine  the  unknown  values  of  y  numerically 
we  have  to  compute  the  coefficients  of  the  matrix 
given  in  (III-19)  at  a  number  of  suitable  control 
points  (xjty  J  on  a ^  equal  to  the  number  of  the 
unknown  p's  .  It  is  known  from  the  potential  theory 
of  double  distributions  that  the  limit  of  the  normal 
derivative  of  the  potential  exists  on  aQ  under  the 
condition  that  the  moment  y  have  a  continuous  second 
derivative  with  respect  to  £  and  n  in  a  neighbor¬ 
hood  of  the  point  (x^tyj)  [we  refer  the  reader  to 
Kellogg,  1929,  pp.  168]  for  more  information.  Accord¬ 
ingly,  the  control  points  should  be  taken  inside  the 
elements  rather  than  at  the  nodal  points  for  our 
analysis. 

In  order  to  be  able  to  compute  the  values  of  y 
near  the  forward  end  of  the  body  accurately,  a  different 
perturbation  scheme  should  be  used  to  solve  the  non¬ 
linear  problem  formulated  in  Chapter  I  near  that  edge. 

We  thereby  obtain  a  local  solution  that  complements 
the  solution  we  have,  and  will  be  matched  with  it 
in  the  next  section  [see  Van  Dyke,  1964,  4.9  for  more 
discussion  concerning  this  method  of  solution] .  Instead, 
we  will  assume  that  the  behavior  of  y  near  to  this 
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edge  can  be  expressed  by  the  following  relation 


Furthermore/  to  simplify  the  numerical  computations 
this  relation  is  approximated  by  a  straight  line, 
which  bounds  the  same  area  as  that  given  by  the  above 
equation,  namely, 

VM+l  -  Vm/3 

Obviously,  there  will  be  a  discontinuity  in  the  pres¬ 
sure  jump  at  this  forward  end  of  the  body.  To  obtain 
the  forces,  the  integration  of  the  pressure  jump  has 
to  be  extended  to  a  point  just  off  this  end. 

The  control  points  (xi,y^)  will  be  taken  as  shown 
in  the  figure  below. 


i 

* 
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NOW/  to  compute  the  coefficients  of  the  matrix 
for  these  control  points,  we  need  some  efficient 
means  of  evaluating  the  double  intrgrals  in  (C-ll) , 
which  we  will  denote  by  Dm,n  (atb)  ,  and  the  inte 
grals  in  (C-15) ,  which  we  denote  by  Im*n  (a,b) 


The  Dm*n  (atb)  Integrals: 


These  integrals  can  be  written  in  the  form, 


Dm,n  (a,b) 


F(a,b; p, \) 


d\ 


We  will  approximate  the  integral  with  respect  to  p 
by  using  the  integration  formula  for  exponential 
integrals  given  by  Abramowitz  (li'68,  25,  4.45).  The 

result  is  the  following: 

n 

Dm'n  (a,b)  *  l  •  Fjfatbsp^)  +  Rn  ,  (D-2) 

iml 

where  p ^  is  the  ith  zero  of  the  Laguerre  polynomial 
bn(p)  ,  w.  is  the  weight  function  given  in  tables 
by  Abramowitz,  and  R n  is  the  truncation  error, 

*n  ’  tMjT  F\ ■  <°  <  «  <  ‘>  ■ 

The  function  F.(tte; p J  is  now  defined  by 
1  t 

00 

Fj(atbiP^)  *  F  (a,b;p^,\)  d\ 

The  general  behavior  of  the  integrand  F(t,8iP^t X) 
is  as  shown  in  the  figure  below. 
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This  integral  is  evaluated  numerically  as  follows: 


F(a,b; p XJ  d\  + 


f1  X2 

*  7T  F(a,b;  p.,y)  d\x  ,  (D-3) 

j  o  u  * 


where  X*  is  the  value  of  X  at  which  the  function 
F(\)  starts  decreasing  monotonically  toward  zero. 
Each  of  the  above  integrals  in  (D-3)  is  computed 
numerically  for  each  value  of  by  using  -he  tra¬ 
pezoidal  rule. 

| 

The  Im,n  (a.bJ  Integrals: 


r  X  i  fX2 

L  i  m 


These  integrals  can  be  written  in  the  form 
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Im,n  (atb)  -  ^  F2(a,bi\)  d\ 


(D— 4 ) 


The  general  behavior  of  the  integrand  F2(a,b;\)  is  as 
shown  in  the  figure  below 


Numerically,  the  improper  upper  limit  of  inte¬ 
gration  will  be  replaced  by  a  finite  value  A0  .  A 
simple  error  expression  E(a,b;\0)  can  be  derived  as 
a  function  of  the  parameters  a,b  and  A0  as  follows: 


E(a,b; \o) 


F2(a,b;\)  d\ 


The  function  F2(atb;\)  in  the  above  integral  can 
be  simplified  considerably  if  we  consider  values  of 
A  where  A  >>  1  .  Consequently  the  above  integral 

can  be  approximated  analytically.  This  expression  can 
be  used  to  find  values  of  Ao  that  achieve  the 
desired  numerical  accuracy.  Knowing  A0  ,  the  inte¬ 
gral  in  (D-4)  can  be  evaluated  by  the  use  of  the 
trapezoidal  rule. 


